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Multiscale phenomena which include several processes oc- 
curing simultaneously at diflterent length scales and exchang- 
ing energy with each other, are widespread in magnetism. 
These phenomena often govern the magnetization reversal 
dynamics, and their correct modeling is important. In the 
present paper, we propose an approach to multiscale mod- 
eling of magnets, applying the ideas of coarse graining. We 
have analyzed the choice of the weighting function used in 
coarse graining, and propose an optimal form for this func- 
tion. Simple tests provide evidence that this approach may 
be useful for modeling of realistic magnetic systems. 
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A large number of phenomena taking place in mag- 
nets include processes occuring simultaneously at differ- 
ent length scales. A good example is the magnetization 
reversal in a macroscopic piece of magnetic material pos- 
sessing different kinds of defects, voids, surfaces etc. The 
reversal starts by a nucleation of a domain with the mag- 
netization opposite to the initial direction. As a rule, the 
nucleation happens near defects, where spins can be frus- 
trated. Here, the different length scales involved can be 
clearly identified. First, there is the microscopic scale 
with a characteristic length of the order of several inter- 
atomic distances (several tens of angstroms) , which corre- 
sponds to the region of spin frustration and contains the 
microstructure in the vicinity of the defect. Next, there is 
a "micromagnetic" length scale (of order of the domain 
wall width, several thousands of angstroms) where the 
formation of the general structure of the nucleus takes 
place. And, finally, the truly macroscopic length scale 
(of order of several microns or even millimeters) , where 
the magnons created in the course of the reversal prop- 
agate. These magnons play an important role in energy 
transfer and sometimes can initiate the magnetization 
reversal in other areas of the sample ||^ . A similar picture 
of several interacting length scales appears in many situ- 
ations, such as the breakthrough of a domain wall pinned 
by a defect influence of the surface on the magnetic 
structure of the core of a magnetic particle Q etc. 

Processes of this type, called multiscale processes, are 
receiving considerable attention nowadays. Along with 
very interesting and rich physics, these are the very pro- 
cesses which govern the switching behavior of magnetic 
systems (coercive field, switching time, etc.), so that an 
adequate understanding of multiscale phenomena is of 
paramount importance for development and creation of 
new magnetic storage media. Micromagnetic simulations 



can provide a realistic description of the processes tak- 
ing place at both micromagnetic and macroscopic length 
scales. On the other hand, microscopic inhomogeneities 
require atomistic simulations (e.g., spin dynamics mod- 
eling 1^ is an adequate tool for materials where spins are 
well localized at the sites of the crystalline lattice). How- 
ever, for the description of real systems, all the length 
scales should be coupled, i.e. modeled simultaneously 
and seamlessly, with the possibility of energy transfer 
between them. 

A simple scheme, where the region of micromagnetic 
simulations is just attached to the region of atomic spin 
simulations, does not give a satisfactory description of the 
energy transfer between the lengthscales. The problem 
is that micromagnetics does not take into account atom- 
istic degrees of freedom; even reduction of the computa- 
tional mesh in micromagnetic modeling down to atomic 
scale does not describe the short-wavelength excitations 
properly, contrary to spin dynamics simulations. The 
artificial sharp boundary appearing between the micro- 
magnetic and atomistic regions leads to the unphysical 
scattering of magnons transferring the energy from one 
region to the other. To couple length scales correctly, 
some transition region between the micromagnetic and 
atomic regions is necessary, which allows for gradual ex- 
clusion of short-wavelength modes until they die out in 
the regions far from the defect. 

Similar difficulties arise in modeling dynamics of struc- 
tural defects in crystals (dislocation motion, crack prop- 
agation, etc.). The coupling of lengthscales for this class 
of problems has been extensively discussed in the liter- 
ature One of the most promising approaches to 
the solution of this problem is coarse-grained molecular 
dynamics (CGMD), developed in Ref. In the present 
paper, we apply the ideas of CGMD to magnetic ma- 
terials modeling. We propose a computational scheme 
which employs basic concepts of nonequilibrium statis- 
tical mechanics to couple micromagnetism and the dy- 
namic modeling of classical spins. We identify the key 
problems arising in the course of implementation of this 
scheme and their possible solutions. 

To explain the idea of the method, let us consider a 
ferromagnet where some magnetic inhomogeneity (de- 
fect, interface, etc.) is present. We describe it as a sys- 
tem of classical magnetic moments of fixed length 
M = |M^|, located at the /i-th site of the crystalline 
lattice (greek indices enumerate the lattice sites). We 
assume that the system is described by the Hamilton 
function 
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n = n° + v, 



(1) 



where Ti!^ describes the isotropic exchange interaction, 
while V Ti.'^ represents ah the other, much weaker in- 
teractions present in the system, and Jf^^, is the exchange 
integral between the sites fi and i'. Let us focus on the 
region far from the defect, where the amplitudes of short- 
wavelength excitations are small (the region where their 
amplitudes are not small should be treated completely 
atomistically). Suppose, a computational grid is defined 
in this region, so that the magnetic state is described 
with required precision by defining the magnetization di- 
rection at the grid nodes: = (MJ,Mj',M|) at the 
j-th node (below, latin indices enumerate the nodes of 
the mesh). These data represent the necessary number 
of large-scale degrees of freedom, and exact modeling of 
all the atomic-scale modes would be too excessive (and 
too expensive resoursewise). However, the atomic-scale 
modes can not be excluded completely, since their cu- 
mulative effect (dissipation, energy transfer, etc.) is not 
negligible, and should be taken into account. 

This problem can be solved by invoking statistical 
physics, i.e. under the assumption of ergodicity the ex- 
act dynamics of short-scale modes can be replaced by 
their statistical distribution. For many relevant situa- 
tions it has been shown |^ that the short-scale modes re- 
lax almost immediately to the state of local quasiequilib- 
rium determined uniquely by large-scale modes, so that 
knowledge of only the large-scale parameters ( "gross vari- 
ables" ) determines the dynamics of the system with nec- 
essary rigor. This description of the system, where only 
large-scale modes are essential (while short-scale modes 
are "slaved" due to the requirement of the local equi- 
librium), is often referred to as the coarse-grained de- 
scription. The theory of local quasiequilibrium states has 
been developed in the 1950s-1960s, and a number of ap- 
proaches exist ^0 11 ijjl^. In the following, we will use 
the convenient modification of the nonequilibrium sta- 
tistical operator (NSO) method 0. For simphcity, we 
do not consider the dissipation processes and quantum 
spin effects, but in principle they can also be included 
following Ref. @. 

It is convenient to introduce new dynamic variables 
ai.p and a2,/j for the magnetic moments M^t according 
to the following relations: 
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These variables can be expressed via conventional polar 
and asimuthal ip^^ angles as 



1 (61^/2) cos (f>f,, 



Msin (6*^/2) sin 0^. 



(3) 



It can be checked that the variables ai^^ and Q!2,p are 
canonically conjugated, so that Hamiltonian formalism 
can be used, and a computational scheme preserving the 
simplectic structure can be employed. 

The conventional way to develop a coarse-grained de- 
scription of the system is to introduce the large-scale vari- 
ables aj as averages of a^: 
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where f^j is an appropriate weighting function, satis- 
fying the normalization condition "^^f^ij — 1, and lo- 
calized near the node j; proper choice of the weighting 
function f^j will be discussed below in more detail. 

Averaging is a natural way to introduce gross vari- 
ables for a ferromagnet, where the local equilibrium is 
governed by the Hamiltonian TC^, so that all the mag- 
netic moments near the node j are almost parallel 
to Mj. Introduced in this way, a(i^2)j constitute col- 
lective variables of the system |l5| changing slowly with 
respect to quickly relaxing short-scale modes, so they can 
be considered as quasi-integrals of motion (for more de- 
tailed discussion see Refs. ^,|l^). Following standard 
procedures of statistical physics, the integrals of motion 
are included into the distribution function via Legendre 
transformations, i.e. we introduce Legendre multipliers 
Fj and Gj corresponding to the node variables aij and 
a2,j, so that the distribution function can be written as 

p = exp (- J2 ml -T.^=ll P^U.j^i,^ (5) 



2ttM ' 

where — 1/ [kBTp), ks is Boltzmann's constant, is 
the spin temperature at the /i-th site, and Q, as can be 
seen, is the normalization factor, analogous to the sta- 
tistical sum of a canonical ensemble in the equilibrium 
case. Note that smooth variations of the temperature 
across the sample can be taken into account (tempera- 
ture is a Legendre multiplier for the integral of energy), 
but to make the consideration simpler, we neglect it, as- 
suming constant temperature, small enough to satisfy the 
condition (3J^^ S> 1. 

The variables Fj and Gj can be considered as parame- 
ters of a local fictitious field acting on the moments M^. 
Values of the parameters Fj and Gj should be chosen in 
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such a way that the resulting equihbrium values a(i^2),/i 
averaged with the weighting function f^j give the re- 
quired values of a(i,2) j- It is easy to see from Eq. ^ that 
the values F^, producing the required values Q!(i,2),j 
are determined from the equations 
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(6) 
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where J- is an analog of the Gibbs' free energy function. 

Now, having all the information at hand, we can use 
averaging to get the equations of motion for the coarse- 
grained variables a(i_2)j- The underlying dynamics of 
the microscopic variables 2),^ is Hamiltonian, i.e. 
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where H is the Hamilton function (|l|). Using the distribu- 
tion function ^ with the values Fj and G° determined 
from Eq. ^ we obtain: 
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where (...) means averaging with the distribution (|^). 
Note that the equations of motion are nonlocal over the 
node indices even if only local exchange interactions are 
present in the system; this important feature is totally 
missing in the micromagnetic description. 

Direct implementation of the scheme presented above 
can be rather expensive computationally. To make the 
problem easier, we can take into account that these com- 
putations are to be performed only inside the relatively 
narrow "belt" between micromagnetic and atomic-scale 
regions (see Fig. |^), where departures of magnetization 
from equilibrium are already small (otherwise, atomic- 
scale simulations should be used). If the z-axis of the 
co-ordinate frame is aligned with the equilibrium direc- 
tion of magnetization the values of ck(i 2),^ are small, and 
the Hamiltonian TiP can be expanded in terms of a(i_2).p 
up to second order: 



^' = \ E E 

u,v — l,2 fi,^ 
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where the indices u,v = 1,2 denote the variables ai and 
a2- In this case, the distribution (||) is Gaussian, so that 
the integral Q and the function T in (^ can be calculated 
exactly [|l6|. As a result, Eq. determining the values of 
generalized torques Fj and Gj can be written as a set of 
linear equations: 
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where y^, is the inverse of the dynamic matrix D^^.m 
in the linearized exchange Hamiltonian (^), and we used 
the vector T^,, (w = 1,2) to denote both F^ and as 
follows: = Fl, and = G^. 

Hence, the problem of dynamic coupling of length 
scales is reduced to two linear problems, Eqs. |l^ and 
^. Note that the rotation of the co-ordinate frame which 
brings the z-axis into coincidence with the equilibrium 
direction of magnetization makes the dynamic matrix in- 
dependent of a^l^2),j^ so that its inverse can be calculated 
once and stored for further references. However, there is 
a subtlety in inverting D^^.m'- this matrix, being deter- 
mined only by the isotropic exchange interactions, has a 
zero eigenvalue corresponding to a shift of all q:(i^2),/j by 
the same value, or, in other words, the dynamic matrix 
has an eigenvector dP = (1,1,. ..,1) corresponding to 
the zero eigenvalue. It reflects the fact that the exchange 
Hamiltonian (^ is invariant with respect to rotation of 
the system as a whole. Thus, when inverting numerically 
the dynamic matrix, a component corresponding to the 
zero eigenvector cf^ should be excluded. 

In the present form, the essense of coarse-graining be- 
comes especially clear. Imagine that, applying some ficti- 
tious torques Fj and Gj we brought the system into such 
a state that the equilibrium values of the large-scale vari- 
ables are a(i^2).j- Then, the atomic magnetic moments 
(i-c, the microscopic variables a(i,2)./^) move in such 
a way that, after the stage of quick relaxation is finished, 
their positions minimize the total energy of the system 
with respect to the constraints imposed by the torques 
Fj, Gj. 

The last but not least problem is an appropriate choice 
of the weighting function Z^.^: it can be verified that 
an arbitrary function does not automatically give mean- 
ingful results. To study this question, let us inspect 
closely the basic idea of the coarse-grained description. 
For a general system consisting of N microscopic mo- 
ments, an exact description of the system's dynamics re- 
quires knowledge of all 2N microscopic canonical vari- 
ables a(i 2).^- However, we expect that under certain 
conditions (which are yet to be formulated), the sys- 
tem can be described with reasonable precision using the 
much smaller set of 2L gross variables a(i^2),j- In other 
words, we expect that under certain approximations, we 
can define such a set of variables a(i,2),j that allows speci- 
fication of all microscopic variables a(i,2),/i with sufficient 
precision. In particular, it can be shown |^ that, for a 
given weighting function /^.j, an optimal (in the least- 
square sense) restoration of microvariables is achieved via 
linear transformation 



E^. 
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where Nj^^ = Y,^f^,k{J2^, f„.jU,k) ^- Thus, the prob- 
lem is to find such a function f^j^j which would make the 
restoration ( pl| ) as accurate as possible. 

In the coarse-grained region, where the linearized 
Hamiltonian (^) can be used, the system's dynamics can 
be represented as an independent motion of different nor- 
mal modes (eigenvectors of the Hamiltonian (H)), and the 
choice of the 2L gross variables becomes obvious: they 
should be amplitudes of the normal modes corresponding 
to lowest L eigenfrequencies. This choice gives an almost 
complete description of the system provided that the am- 
plitudes of all the other modes, which are excluded from 
consideration, are much smaller. The frequencies of the 
excluded modes are much larger, so that their dynam- 
ics is much faster, and they can relax to local equilib- 
rium quickly in comparison with the adiabatically slowly 
varying gross variables. 

This choice of the gross variables is in correspondence 
with the dynamical approach to nonequilibrium sta- 
tistical mechanics. When considering motion of the sys- 
tem subjected to some small rapidly varying perturbation 
(provided, in our case, by the fast excluded modes), the 
well-known basic problem is to exclude from the solu- 
tion so-called secular terms, which appear due to entan- 
glement of slow and fast motions in the system. After 
slow and fast modes are properly separated, the stan- 
dard procedure of averaging can be performed. The use 
of the lowest-frequency normal modes as gross variables 
allows elimination of the secular terms in the coarse- 
grained equations of motion. An analogous procedure 
can be identifycd also in the approaches of Zwanzig ||l^] 
and Mori 

Thus, one possible recipe is to use the weighting func- 
tion ffj_j = d^, where is one of the lowest-frequency 
eigenvectors of the dynamic matrix. In real calculations 
it could be inconvenient to work with delocalized collec- 
tive modes. It can be shown, that an equally accurate 
coarse-graining can be achieved with any orthonormal 
combination of the eigenvectors , i.e. f^j — J^k ^l^'^i 
is an equally good weighting function if 

T.cici=s,,r T.^lCi,^Sk,k'. (12) 

k j 

The function /^j- can be made well-localized using an 
appropriate set of C^. 

As a specific example, we performed coarse-grained 
modeling of magnons in a 1-D ferromagnetic chain con- 
sisting of N classical magnetic moments with nearest- 
neighbor and next-nearest-neighbor exchange interac- 
tions. The corresponding Hamilton function is 

?i = ^J°M^(M^+i+M^_i) (13) 
+7j;]M^(M,,+2+M^_2), 



and periodic boundary conditions are used. A computa- 
tional mesh is imposed, consisting of L nodes, and two 
gross variables aij and are attributed to each of L 
nodes. These two gross variables provide a coarse-grained 
description of a block containing n — N / L individual mo- 
ments. 

Several weighting functions have been used. One pos- 
sible choice, giving an exact magnon spectrum, is 

(0) ^ 1 sin7r(^- jn)/n 
■^^'J TV tan7r(^- jn)/iV' ' 

which corresponds to an orthonormal combination 
of exact normal modes with coefhcients Cj, — 
(1/a/]V) exp (— ij'/cn) where i = T (so it is a discrete 
analog of the Nyquist-Shannon uniform sampling). How- 
ever, this function, due to its long tails, is not convenient 
for computations, and we tested its rescaled cutoff mod- 
ification: 

f^l^Aflv Ai-J«<«, (15) 
4^] =0, ^ - jn > n, 

where A is the normalization factor necessary to satisfy 
the condition /^j- = 1. A second form of the weight- 
ing function 

4^] = 1/"' ~ j'^ ^ (16) 
fl,j =0, II- ]n> n, 

although very far from optimal, can be used for crude 
semi-qualitative computations, so we also tested its per- 
formance. 

The results of our tests, the dispersion curves uj{k) 
for magnons with different wave vectors k, and their 
group velocities c(k) — dio/dk, are calculated with dif- 
ferent weighting functions, as shown in Fig. ^j. Disper- 
sion curves uj{k) describe the magnon dynamics in the 
chain, while the group velocity curve c{k) characterizes 
the propagation of magnons (the latter should be tested 
separately since a good approximation for LLj{k) does not 
necessarily imply a good approximation of its deriva- 
tive). For comparison, exact curves are presented, along 
with the results of micromagnetic calculations. The data 
points for fc = and k = 2tt /n are not shown because the 
dynamic matrix formally has a singularity at these val- 
ues of wave vector. For very long- wavelength magnons all 
types of modeling work rather well, but for shorter wave 
lengths the differences are large. The coarse-grained de- 
scription is better than the micromagnetic even for the 

(2) 

worst function ■ . For the appropriately chosen weight- 
ing function /^^j , in spite of the cutoff, the difference with 
the exact solution is very small, even at maximal allowed 
wave vector values. 

Summarizing, we propose an approach to modeling of 
multiscale processes in magnets, applying the ideas of 
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coarse grained molecular dynamics ||^ to magnetic mod- 
eling. The scheme proposed employs basic concepts of 
nonequilibrium statistical mechanics to couple length- 
scales. We discussed possible implementation of this ap- 
proach, paying particular attention to the problem of the 
correct choice of the weighting function used in coarse 
graining. Simple tests verify our conclusions and evi- 
dence that this approach can be applicable to larger and 
more complicated systems. 
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FIG. 1. In the vicinity of the defect (the defect is depicted 
as a dashed circle in the center) , individual magnetic moments 
can depart significantly from ideal ferromagnetic order, be- 
cause of lack of co-ordination, and due to changes in exchange 
and anisotropy constants. This region can be described using 
atomistic modeling of individual moments. In the region far 
from the defect, spins are almost parallel to each other, and 
a micromagnetic description is appropriate. To couple these 
regions correctly, a transition region is included, where the 
coarse grained description is used. 

FIG. 2. Calculated dispersion curves uj{k) (above) for 
magnons with different wave vectors k, and their group ve- 
locities c(fc) (below). Wave vector values k are shown in units 
of fco — 27r/n. Calculations have been done for the chain con- 
sisting of A*' = 1024 individual moments with the mesh con- 
taining L — 32 nodes. The value of the next-nearest exchange 
integral is J = 1 and the ratio 7 = 0.2 of nearest-neighbor 
to next-nearest-neighbor exchanges are used. The results for 
micromagnetic modeling (MM), and for two versions of coarse 
graining with different weighting functions (CG 1 for fj^j and 

CG 2 for /^^] ) are presented. Exact curves are shown for com- 
parison. 
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